On the phase behavior of hard aspherical particles 
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We use numerical simulations to understand how random deviations from the ideal spherical 
shape affect the ability of hard particles to form fee crystalline structures. Using a system of hard 
spheres as a reference, we determine the fluid-solid coexistence pressures of both shape-polydisperse 
and monodisperse systems of aspherical hard particles. We find that when particles are sufficiently 
isotropic, the coexistence pressure can be predicted from a linear relation involving the product of 
two simple geometric parameters characterizing the asphericity of the particles. Finally, our results 
allow us to gain direct insight into the crystallizability limits of these systems by rationalizing 
empirical data obtained for analogous monodisperse systems 



I. INTRODUCTION 



Understanding how objects pack or can be designed to 
tile the three dimensional space is a fundamental opti- 
mization problem that has important practical applica- 
tions that range from the macroscopic, such as the effi- 
cient storage of grains, to the microscopic: the fabrication 
of band-gap photonic materials. Mathematicians have 
been intrigued by such problems for centuries; namely 
since Kepler's 1611 essay On the Six-cornered Snowflakes, 
but it wasn't until recently, with the advent of nanotech- 
nology and the explosion of molecular biology, that prob- 
lems of packing and self-assembly of nanoscopic compo- 
nents gained tremendous traction in the broader scientific 
community. Recent advances in synthesis of nanopar- 
ticles |2h10| allowed for unprecedented control over the 
shape and surface chemistry of colloidal particles, thus 
providing an unlimited number of building blocks whose 
spontaneous aggregation could lead to the formation of 
an unprecedented variety of structures with potentially 
novel functional, mechanical, and optical properties. Un- 
like most of the work on particle crystallization and self- 
assembly that in the last decade has focused on monodis- 
perse [ll| or polydisperse [H, EH systems of spherical or 
regularly-shaped particles (see also 0, [Ur{24[ and ref- 
erences therein), surprisingly little or nothing has been 
done theoretically to understand the packing of irregu- 
larly shaped particles. Indeed, there are several impor- 
tant cases in which the shape of the single components 
cannot be tailored at will, yet, an efficient packing, or an 
understanding of the physical properties of these densely 
compressed systems, is highly desirable. An example of 
an outstanding problem in this category is the storage of 
grains (25| . 

Here we seek to gain insight into the crystallization of 
hard aspherical particles. We want to understand how 
particle geometric features can be related to their abil- 
ity to orderly pack into three dimensional periodic struc- 
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tures, and especially to identify under what conditions 
they cease to do so. We have recently reported [l[ that it 
is possible to empirically relate particle geometry to crys- 
tallizability (intended as the tendency of a component to 
crystallize) by using two simple geometric parameters. 
The first is the particle asphericity A, defined in terms 
of the surface to volume ratio of a particle a p = A p /V p 
with respect to that of a sphere of diameter a, a s = 6/cr, 
as A = 1 — a s /a p . The second parameter, q, is related to 
the orientational symmetry of the particle. It is used to 
describe the asphericity of random walks [26[ , and it is 
obtained by combining invariants of the particle inertia 

tensor I tJ as q = Ei<j'( A i ~ V ' where Ai ' 

with i = 1, 2, 3, is an eigenvalue of Uj. 

In this paper we try to rationalize those empirical re- 
sults by computing how random perturbations from the 
ideal spherical shape affect the fluid-solid coexistence 
pressure of monodisperse and shapc-polydispcrse systems 
of hard aspherical particles. 



II. METHOD 

In order to generate a statistical ensemble of hard as- 
pherical particles, we developed a simple model [l[ that 
guarantees a certain degree of control over the particle 
shape. Each particle is built by setting the center of At 
(4 < Ah < 12) spheres of diameter a at random posi- 
tions on the surface a sphere of diameter oo < a, allowed 
to overlap freely. Polydisperse systems of such aspherical 
particles were generated by choosing specific values of Af, 
and (To and allowing each particle in the system to arise 
from a different random collection of sphere positions. 
Deviations from the spherical shape can be conveniently 
controlled by varying o~q and A&. 

All dimensions of the resulting particles are then scaled 
so that the volume V p of each particle is equal to that of 
a sphere with radius a, V s = \ko- 2 . The resulting sphere 
diameters after this rescaling are denoted o~r. For <7q = 
one recovers the spherical limit, and as Co increases, 
particles develop larger and larger shape distortions. In 
a similar fashion, large values of A;, result in a bumpy but 
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overall isotropic particle, whereas small values of Nf, tend 
to generate very anisotropic shapes. Any two particles i 
and j interact via a hard repulsive potential defined as 



Ui 



if 



r t \> o-r Vs G i , Vt G j 



oo otherwise 



(1) 



where s and t run over all spheres of rescaled diameter an 
constituting particle i and particle j respectively. Exper- 
imental realizations of colloidal particles similar to ours 
could be generated using the approach described in ref- 
erence [7H9| to create nonspherical particles with tunable 
shapes. 

In order to determine the fluid-solid coexistence pres- 
sures for both polydisperse and monodispcrsc systems of 
asphcrical particle, we used the method of direct fluid- 
solid coexistence simulation described in 27 1. 1024 hard 
particles were placed in an FCC crystal lattice (at volume 
density As — 0.545, the hard sphere crystal coexistence 
density [28[ ) centered in a box of d imensions L x x Ly x L z . 
An FCC crystal was chosen based upon our previous 
work, which indicated that when monodispcrsc asphcr- 
ical systems formed crystals, they tended to be (apart 
from a few exceptions) FCC; we have found no evidence 
supporting the use of any other crystal geometry. The 
dimensions of the box were chosen such that L x and L y 
were just large enough to accomodate the FCC crystal, 
and L z was roughly four times larger. The crystal lat- 
tice was chosen such that the extension of the crystal 
in the z-direction was roughly twice that in the x- and 
y-dircctions, in order to increase the separation between 
the two fluid-solid interfaces in the system. This crystal 
lattice was placed into equilibrium with a fluid of 1024 
particles at hard sphere fluid coexistence volume density 
pi ~ 0.495. The fluid and crystal were both briefly al- 
lowed to relax in order to relieve any overlap introduced 
by the "bumpiness" of the particles and to allow the fluid 
to come fully into contact with the solid interface. 

Monte Carlo simulations were then run in the NPT 
ensemble, where the three box dimensions L Xl L y , and 
L z were allowed to fluctuate independently. 

The number of crystalline particles Nx in the system 
was monitored as a function of Monte Carlo step; this 
quantity was determined using the standard spherical- 
harmonics based bond order parameter g6 [l^jHoj. Direct 
simulation is known to have a few caveats: slow equili- 
bration, non-insignificant finite size effects, dependence 
of the surface free en ergy on the specific face the crystal 
exposes to the fluid [27|. Nevertheless, for this specific 
system, we find this direct method to be more reliable 
than the two-step thermodynamic integration scheme de- 
scribed in 
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Our empirical results on the crystallizability of systems 
of monodisperse hard aspherical particles [l|, obtained 
by slowly compressing an ensemble of 487 different shape 
realizations (see [l[ for details) are summarized in Fig. Q] 
and indicate the existence of a clear boundary between 
particles that crystallize and particles that do not. A 




FIG. 1: Crystallizability of monodisperse hard aspherical par- 
ticles characterized in terms of two shape parameters A and q. 
Filled circles indicate particles that easily crystallized, while 
open squares indicate particles that did not. The solid line is a 
hand-constructed fit to the crystallizability limit; the dashed 
line represents a prediction of the crystallizability limit from 
the current work; see below. Figure adapted from 0|. 



roughly inverse relationship is clearly evident; particles 
with large A must have very small q in order to have a 
hope of crystallization, and vice-versa. This result pro- 
vides a very useful way of predicting whether a particu- 
lar particle shape can pack into a crystalline structure by 
simply measuring the experimentally accessible A and q. 
As this model is intended to describe randomly shaped 
particles, the diagram docs not include the results for 
particles designed with very specific shapes such as rods, 
plates or regular polyhedric geometries that are known to 
crystallize. These particular cases would generate sharp 
peaks around specific values of q, and are purposefully ex- 
cluded from this study. The solid line is a guide to the eye 
and has the functional form A(q) = 0.023+ l/(170g- 10). 



III. POLYDISPERSE SYSTEMS 

The question remains as to why the systems behave 
this way; this is the focus of the current work. To get 
a better understanding of the crystallizability limit, we 
compute the coexistence pressure of systems of shapc- 
polydisperse hard particles as a function of the geomet- 
ric parameters A and q via the direct fluid-solid sim- 
ulation scheme described above. Incorporating shape- 
polydispcrsity in these systems is however not a trivial 
matter and requires some discussion. Unlike the common 
notion of polydispcrsity of spherical particles for which 
any particle size can be indifferently used as a reference, 
each particle shape whose values of A and q lead to a 
point below the crystallizability boundary shown in Fig.Q] 
could be used as a reference shape for this study. The 
issue is that the phase behavior could be very much de- 
pendent on the specific choice of A and q, making the 
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FIG. 2: Coexistence pressure vs. Aq for shape-polydisperse 
systems of hard aspherical particles. The solid line is a fit 
to the data for Nb > 6. Inset: Distribution of q values for 
polydisperse systems with Nb = 4, 6, 8, and 12. Recall that q 
is independent of cto- The distribution is significantly larger 
and broader for Nb = 4 than for the larger values of Nb- 

problem intractable. We therefore adopted a pragmatic 
approach to describing shape-polydispersity that has a 
natural experimental counterpart Q. Namely, we con- 
sider the distributions of qs and As arising when con- 
structing our particle using a given number of spheres 
Nb and a given value of <7o- Particles with large val- 
ues of Nb generate narrow distributions shifted towards 
small values of q, while small values of Nb (Nb > 4) re- 
sult in wide distributions peaked over large values of q 
(see Figure [2Ja)). When cr one recovers the hard 
sphere limit for any value of Nb, and as <jq increases, 
larger and larger values of A will be sampled (without 
affecting the q distributions) until approaching the crys- 
tallizability boundary. 

We systematically analyzed hard particles constructed 
with Nb =4, 5, 6, 8, and 12, and values of <tq 6 
[0.05,0.30]. 

All data points have been computed using two differ- 
ent random sets of particle geometries from the given 
distribution. The difference between the results from 
the two independent sets is smaller than the error bar 
associated with the numerical scheme. Figure [5] shows 
our results. Remarkably, we find that the coexistence 
pressures, p*(A,q), between fluid and solid for values 
of Nb > 6 can be made to collapse into a simple lin- 
ear dependence on Aq which can be fitted to p*(A, q) = 
99lAq + 12; note that the y-intercept is near the hard 
sphere coexistence pressure of pj^s — H-? EH 7 as it 
should be. Our results thus suggest that p*(A,q) can 
be expressed in terms of a Taylor expansion in the quan- 
tity Aq: p (. A 'i)-Pns _ a A„ -)-..., Note, however, that 

Pus 

significant but systematic deviations from the linear fit 
are visible for systems of particles obtained with Nb = 4 
and Nb = 5. These values of Nb correspond to parti- 
cles with relatively large values of q and rather broad 
distributions (see Figure [2] inset), suggesting that higher 
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FIG. 3: Coexistence pressure vs. Aq for monodisperse sys- 
tems of hard aspherical particles with Nb — 6, ao — 0.2. The 
line is a fit obtained from the data presented in Figure [2] for 
polydisperse systems with Nb = 6, 8, and 12. The cross in- 
dicates the position of the polydisperse system with Nb = 6, 
00 = 0.2. Points at the extremes of the q distribution are 
shaded with a lighter color. 

order terms in the expansion may be necessary. Further- 
more, for such wide distributions fractionation between 
isotropic and anisotropic particles may become an im- 
portant factor. 



IV. APPLICATION TO MONODISPERSE 
SYSTEMS 

The existence of a master curve for the coexistence 
pressure of significantly different particle distributions 
(provided Nb > 6), both in terms of their average value 
and their width, is quite remarkable, and provides a very 
compact and elegant way of expressing the equilibrium 
properties of aspherical systems, thus significantly reduc- 
ing the number of variables required to describe these 
complex systems. The results from the above study are 
very robust. In fact, we next considered 10 monodisperse 
systems generated by a single set of parameters: Nb = 6, 
(To = 0.2. Each of these systems is characterized, by def- 
inition, by an infinitely sharp P(q), and individual par- 
ticle geometries were selected to mostly cover the range 
of Aq values explored above, while keeping A roughly 
constant because for a given Nb and a$, the distribution 
of A is quite narrow. A plot of these data, shown in 
Figure [3J indicates that the relationship adapted unal- 
tered from the polydisperse case is a good predictor of 
coexistence pressure for the monodisperse case as well. 
This is a clear indication that, as long as fractionation 
effects do not come into play, the master curve is overall 
independent of distribution width. This behavior can be 
readily linked to the plastic nature of the crystals formed 
by aspherical monodisperse systems [l[ for comparable 
values of q and A, i.e. each particle can freely rotate 
around their lattice sites, and can be understood quali- 
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tatively by realizing that any two identical neighboring 
particles will typically be misoriented and face each other 
with a random region of their respective surfaces. As 
a consequence, their mutual interaction becomes hardly 
distinguishable, on average, from that of two particles 
with different shape, thus explaining why monodisperse 
and polydispersc systems may indeed present analogous 
equilibrium properties. Note, however, the deviations 
from the predicted pattern at the extreme ends of Aq, 
corresponding to particle shapes that have values of q 
significantly different from the average q of the distribu- 
tion. The deviation at large q is expected as the analo- 
gous behavior is observed for polydisperse systems; as 
particles become more anisotropic their rotational de- 
grees of freedom are reduced until they perfectly align 
in the q — > 1 (rod) limit, making the averaged random- 
shape argument described above inappropriate. A bit 
more surprising is the deviation for very small values of 
q, for which we find that the polydisperse system with 
A r 6 = 12 nicely follows the master curve. To rationalize 
this behavior one has to realize that for small values of 
TVb, when q — > 0, particles tend to acquire rather symmet- 
ric and specific geometries, such as platonic solids, which 
also require oricntational ordering to tile the space as 
soon as A becomes sufficiently large. These specific par- 
ticle shapes dominate the shape space for q ~ and small 
values of Nb, and, as discussed in lead to deviations 
from the inverse power law behavior. 

Notice that no coexistence pressures above p* (A, q) ~ 
21 are reported. This is interesting because this value 
is close to the fluid pressure of a system of hard spheres 
at the glass transition density [32], Pg = 22.6, and it 
is reasonable to assume that pa sets an upper bound 
to the largest accessible fluid coexistence pressure for 
hard aspherical particles obtained with direct fluid/solid 
sampling. In fact, this method is clearly susceptible to 
anomalies in system kinetics, thus making coexistence 
measurements at pressures larger than the one here re- 
ported quite cumbersome and somewhat unreliable. Nev- 
ertheless, this provides a way of rationalizing the empir- 
ical data, because if we assume that no aspherical hard 



particle will easily crystallize above pa, we can readily 
link the pressure master curve p* (A, q) to the inverse 
power-law behavior obtained for monodisperse aspheri- 
cal systems. By plugging pa into p*(A, q), we obtain the 
crystallizability boundary Aq ~ 0.011; this equation is 
plotted on Figure [T] along with the original line presented 
in 0. 

Although it is not a perfect division, likely due to many 
factors (including finite-size effects), this simple predic- 
tion provides a good dividing line which would predict 
the packing into fee crystals of about 80% of all systems 
tested. Given the uncertainties associated to the data 
points in [l[ , we believe this is a remarkably good rate of 
success, and provides an elegant physical reason for the 
inverse power law relationship found therein. 

V. CONCLUSIONS 

To conclude, our data provide a simple way of ratio- 
nalizing the crystallizability limit observed for systems of 
monodisperse hard aspherical particles and generalize the 
results to include systems containing non-identical hard 
aspherical particles. Our results show that the coexis- 
tence pressure of systems of monodisperse and shape- 
polydisperse hard particles with fee crystals is remark- 
ably well-described, within the limits described in the 
text, by a simple linear relation in Aq, p ~ PffS = aAq, 

Phs 

and indicate precise limits for the manufacture of nano- 
componcnts expected to crystallize. Our study is by no 
means exhaustive of all possible shapes, but it represents 
a first major step towards a statistical understanding of 
the role of shape disorder in particle organization. Our 
theoretical predictions should be readily testable experi- 
mentally. 
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